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ABSTRACT 

We present the serendipitous discovery of a young stellar cluster in the Galactic disk at 1=12°. 
Using Keck/NIRSPEC, we obtained high- and low-resolution spectroscopy of several stars in the 
cluster, and we identified one red supergiant and two blue supergiants. The radial velocity of the red 
supergiant provides a kinematic cluster distance of 4.7±0.4 kpc, implying luminosities of the stars 
consistent with their spectral types. Together with the known Wolf-Rayet star located 2.'4 from the 
cluster center, the presence of the red supergiant and the blue supergiants suggests a cluster age of 
6 — 8 Myr, and an initial mass of > 2000 Mq. Several stars in the cluster are coincident with X-ray 
sources, including the blue supergiants and the Wolf-Rayet star. This is indicative of a high binary 
fraction, and is reminiscent of the massive young cluster Westerlund 1. The cluster is coincident 
with two supernova remnants, SNR G12.72— 0.0 and G12.82— 0.02, and the highly magnetized pulsar 
associated with the TeV 7-ray source HESS J1813— 178. The mixture of spectral types suggests that 
the progenitors of these objects had initial masses of 20 — 30 Mq. 
Subject headings: stars: evolution — infrared: stars — X-rays: stars — stars: supernovae: general 

1. INTRODUCTION 



Young stellar clusters are important tools for investi- 
gating the structure, the chemical enrichment, and the 
current star formation of a galaxy. Furthermore, they 
are natural laboratories to study the evolution of massive 
stars. Massive stars explode as supernovae creating neu- 
tron stars and black-holes, and are believed to be sources 
of 7-ray bursts, which are the most energetic explosions 
in the Universe. Because of their short lifespans, mas- 
sive stars are rare and predominantly observed in young 
(few Myrs) massive (~ 10** Mq) clusters. Only a hand- 
ful of massive stellar clusters are known in the Galaxy: 
important examples inclu de Arches, Q uintuplet, Wester- 
lund 1, and Westerlund 2 (|FigeiJl2008[ ). As these clusters 
lie in the plane, their detection is hampered both by in- 
terstellar extinction and by the relatively low numbers 
of member stars. Indeed, the census is severely incom- 
pletely, as demonstrated by the recent disco very of two 
other massive clusters (R SGCl and RSGC2. lFiger et all 
[200I iDavies et al.l[2007l) . 

We have serendipitously discovered a new stellar 
cluster at RA= 18hl3'"24!31 and DEC=-17°53'30'.'83 
(J2000). Previously undetected, this cluster clearly ap- 
pears as a stellar over-density at infrared wavelengths 
(Fig. [T]), and is also visible at optical wavelength. It 
is located adjacent to the star forming region W33 in 
an extremely rich and complex region. Two super- 
novae remnaiits (SN Rs), G12.82-0.02 and G12. 72-0.0 
(iBrogan e t all 120061): a -7-rav source, HESS J1813-178 
( Aharoni an et al.l |2006[ ) ; and numerous X-ray sources 
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Fig. 1. — i^s-band image of the cluster from 2M ASS, with over- 
layed contours of 90 cm data llWhite et al.l l2005'). Contour levels, 
in mjy beam-\ are: 40, 60, 80, 100, 120, 140, 160, 200. The beam 
size is 24" X 18", FWHM, and the position angle of the major axis 
is along the North-South direction. Stellar identification numbers 
are from Table 1. 
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(|Funk et all l2007t iHelfand et all 120071) have been de- 
tected in its surrounding area. 

In this report, we present near-infrared spectra of four 
cluster members, of which three are identified as early 
B-type stars, and one as a red supergiant (RSG). Also, 
we analyze existing photometric infrared data and X- 
ray data, and find that several bright cluster members 
are associated with an X-ray so urce, one of which i s 
a known Wolf-Rayet (WR) star (|Hadfield et al.ll2007l ). 
From these results, we estimate the age and distance 
of the cluster, which appears associated with the star 
forming region W33 and the high-energy source HESS 
J1813-178. 
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2. OBSERVATIONS AND DATA 

Spectroscopic observations of several candidate clus- 
ter members (see Table 1) were performed at the 
Keck Observatory o n 2008 April 19 using NIRSPEC 
(jMcLean et al.lll998D under program U050NS (P.I. R.M. 
Rich). A spectrum of star #1 {Ks= 3.7 mag) was ob- 
tained with the NIRSPEC-7 filter and the 0"57 x 24" 
slit, covering 1.99 — 2.39 /im at a resolution of R=17000 
with two nodded exposures of IDs each. Low-resolution 
spectra of sources #2 and #3 were taken with the K fil- 
ter and a 0'.'57 slit width, covering 1.9 — 2.35 /.tm at a 
resolution of R=1700. Data re duction was performed as 
described in iFiger et al.l (|2003D . We subtracted pairs of 
nodded frames and flat-fielded them. The distorted 2- 
D spectral traces were rectified onto a linear grid, using 
arc and etalon frames for wavelength calibration. Atmo- 
spheric absorption and instrumental response were re- 
moved by dividing each extracted target spectrum by 
the spectrum of a BIV telluric standard (HD164581). 

Photometric measurements of stellar point sources cov- 
ering the cluster region were obtained from the Vizier 
database. We cross-correlated 2 MASS near-infrared 
measurements (ICutri et all I2003D. Spitz er/GLIMPSE 
mid-infrared data ( Benjamin et al.l |2003( ). as well as 
optical data from th e astrometric catalog NOMAD 
(jZacharias et al.ll2004f) . 

The regio n surrounding the clu ster was observed with 
Chandra bv iHelfand et all (pOOT), who detected 75 X- 
ray sources. We cross-identified the Chandra sources 
with the 2MASS and GLIMPSE catalog using a radius 
of 1'.'5, and found 44 matches. 

3. SPECTRAL AND PHOTOMETRIC ANALYSIS 

In order to confirm that the observed over-density is an 
actual stellar cluster, we analyzed the photometric prop- 
erties of stellar point sources from 2MASS and GLIMPSE 
surveys. As an example, a {J— Kg) vs K^ diagram of 
2MASS point sources within 3f5 from the peak of over- 
density is shown in Fig. [H A well defined cluster se- 
quence appears in the infrared color-magnitude diagram 
(CMD) at about J-Ks^ 1.5 mag and Ks= 4-12 mag. 
Its red color and bright magnitude distinguish the cluster 
sequence from a foreground component. An interstellar 
extinction of ^Ks= 0.83±0.2 mag (Av= 9.1 mag) is esti- 
mated by matching the colors of the observed cluster se- 
quence with a theoretical isochrone (6.3 Myr, solar com- 
posit ion) from the Geneva group (jLeieune fc Schaere^ 
l2001h . and by assuming a powe r-law extinction curve 
A\ oc A~^-^ (jMessineo et al.|[2"005[ ). Because the main se- 
quence is almost a vertical sequence when viewed in the 
near-infrared, this estimate is independent of age. 

A variety of massive stellar objects was detected to- 
ward the new stellar cluster: a RSG star, a WR star, 
several BSGs (Table 1). Star #4 is a.ssocia ted with a 
known WR star (#8 in iHadfield et alll2007h . Star #1, 
a star with Ks= 3.7 mag, 3.5 mag brighter than other 
cluster members, dominates the infrared cluster surface 
brightness. Its high-resolution spectrum presents CO- 
bands in absorption as shown in Fig. [3l A radial ve- 
locity Vlsr — +62 ± 4 kms~^ was measured by cross- 
correlating its spectrum with that of Arcturus after re- 
binning the latter at the same resolution. A determina- 
tion of the spectral type was obtained by comparing the 
CO equivalent width with that of other template stars. 
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Fig. 2.— 2MASS J-K^ vs. K^ CMDs. In the left panel, we show 
the CMD of point sources that are within 3.'5 from the center. In 
the right panel, field stars are plotted, which were taken from an 
annulus of equal area at a r adius of 5'. The vertical li ne indicates 
a solar isochrone of 6.3 Myr IILeieune &: Schaererll200ll) . which was 
shifted to a reddening of j4Ka=0.83 mag and a distance of 4.7 
kpc. The starred symbol indicates the location of the WR star, 
triangles the Chandra point sources, the large dot the RSG star, 
and the circles the BSGs. 
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Fig. 3. — Top panel: high-resolution spectrum of the RSG star in 
the region of the CO band— head feature at 2.293 fini. As a compar- 
ison, a spectrum of a K3 I star (Davics et al. 2008) is also shown 
(dashed line). Bottom panel: low-resolution spect ra of #2, #3, 
and # 11. Helium and hydrogen lines are taken from lNaiarro et al.l 
lITOQl) . 



as explained in iDavies et al.l ()2007t ) . The resulting CO 
equivalent width is consistent with that of a K2 I or an 
M4 III star. From the low-resolution frames, three spec- 
tra were extracted: the two spectra of sources #2 and 
#3, and a third spectrum of source #11, which is 6" 
away from #2, and fell on the slit. The detection of 
Hel lines, at 2.058 /zm and 2.112 ^m, and HI line at 
2.166 /im, and the lack of HcII indicate that #2, #3, 
and # 11 are early B-type (B0-B2) stars (jHanson et al.l 
I1996D . The emission line at 2.058 fim suggests that #2 
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List of 


CANDIDATE MASSIVE STARS. 














ID 


Ra 


Dec 


B 


V 


R 


J 


H 


Kg 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


IDx 


Lx 


Sp 


01 


18:13:22.255 


-17:54:15.55 


14.74 


13.68 


12.73 


5.46 


4.23 


3.79 


5.07 


4.05 


3.57 








K2I 


02 


18:13:23.256 


-17:53:26.54 


20.08 


15.03 


14.50 


8.64 


7.72 


7.22 


6.94 


6.78 


6.63 


6.68 


39 


4.1 


B0.5-2I 


03 


18:13:24.431 


-17:52:56.75 


20.42 


15.74 


14.68 


9.20 


8.25 


7.79 


7.47 


7.38 


7.30 


7.30 


43 


4.8 


B0.5-2I 


04 


18:13:14.188 


-17:53:43.60 


20.48 


17.57 


15.20 


9.62 


8.60 


7.94 


7.34 


6.93 


6.70 


6.35 


24 


80.5 


WN7b 


05 


18:13:23.711 


-17:50:40.41 


19.51 


16.41 


14.94 


9.96 


9.06 


8.56 


8.23 


8.10 


7.99 


8.06 


41 


63.4 




06 


18:13:33.817 


-17:51:48.89 






16.12 


10.29 


9.15 


8.57 


8.00 


7.90 


7.75 


7.75 


58 


4.0 




07 


18:13:22.519 


-17:53:50.14 


18.63 


15.60 


16.57 


10.30 


9.27 


8.66 


8.03 


7.70 


7.49 


7.30 


37 


10.3 




08 


18:13:14.494 


-17:54:01.77 


19.91 


15.47 


15.71 


9.86 


9.12 


8.75 


8.43 


8.33 


8.31 


8.31 


27 


3.8 




09 


18:13:21.911 


-17:51:35.91 


21.17 




16.24 


10.74 


9.80 


9.34 


8.95 


8.90 


8.83 


8.91 


36 


2.8 




10 


18:13:47.561 


-17:57:01.43 


15.32 




13.64 


10.76 


9.89 


9.60 


9.31 


9.34 


9.22 


9.20 


71 


21.0 




11 


18:13:23.620 


-17:53:24.50 








11.74 


10.85 


9.59 


10.16 


10.32 


10.02 








BO.5-2 


Note. — For each star, number 


designat 


ions and coordinates (J2000) 


ire followed by 


magnitudes measured 


in different bands. J,H, 


and Kg measurements are from 2MASS, while the 


magnitudes at 3.6 um, 


4.5 um, 


5.8 um 


and 8 


/xm are 


from GLIMPSE. B, 


V, and R 


associations are taken from the astr 


ometric 
es of Lx 


catalog NOMAD. For sources with X-ray associations 
are obtained for a distance of 4.7 kpc, by assuming 


, number designations (IDx] 
N{H) = 1.6 X 10^2 cm-2, a 


are taken 


fromlHelfand et al.| ^20071). Estimat 


power law 



model, and a photon index of 1.5 l|Townslev et al.ll2006I V X-ray luminosities Lx are given in lO^'^erg s ^. 



is a blue supergiant (BSG). 

4. MASSIVE STARS AND DISTANCE 

From the radial velocity of star #1, we derive a kine- 
matic heliocentric distance of 4.7 ±0-4 kpc by using 
the rotation curve of Brand fc Blitj (Il993il and a solar 
Galactocentric distance of 7.6 kpc (jKothes fc Dougherty! 
I2007D . Since the far distance (10.3 kpc) would make the 
stars over-luminous, we only consider the near distance. 
A spectro-photometric distance of 0.9 kpc or 4.4 kpc is 
inferred when assuming an JVI4 giant or an early K RSG, 
respecti vely, together with the absolute infrared magni- 
tudes of IWainscoat et ahl ()1992D . An interstellar extinc- 
tion ^K^ of 0.48 mag is derived. We conclude that this 
bright star is a K2I RSG star. Star #1 is 0f9 from the 
center. 

By assuming a distance of 4.7 kpc, as calculated for 
the RSG star, an intrin sic H—Kg from 0.11 to 0.27 mag 
(ICrowther et al.l l2006al). and using the extinction law 
bv iMessineo et aLl ( 20051 ). we derive for the WN7b star 
an Ak, from 0.83 to 0.59 mag and an absolute mag- 
nitude Mk from —6.0 to —6.2. This value of Mk is 
consiste nt with those of simil ar WN7 stars in Wester- 
lund 1 (jCrowther et al.ll2006aD . The 2MASS colors and 
magnitudes are consistent with being a cluster member. 
A bolometric correction of —3.5 (Crowthcr et al. 2006a) 
yields a bolometric magnitude of —9.6. Star #4 is 2.'4 
from the center. 

The 2MASS photometric measurements of ^2 and ^3 
support either early nearby B dwarfs or early B super- 
giants at a distance of 3.7 ± 1.7 kpc. We assume the 
ab solute magnitudes and bolo metric corrections given 
bv ICrowther et all (|2006bD and lBibbv et"an (l2008f). and 
cons ider both extinction law s bv HVlessineo et al.l ()2005D 
and llndebetouw et all (|2005f ). St ar #2 has ^k = 0.8 7 
mag, and #3 has ^k5= 0.81 mag (jMessineo et al.|[2"005[ ). 
However, they have the same interstellar extinction as 
the WR star; they are located in the same region of the 
CMDs, and most likely they are BO-2 supergiants, as 
well as members of the cluster. They are within 0.'6 from 
the center. Star #11 has a poor 2MASS photometry; 
with Ks~ 9.59 mag and an extinction A^^— 1.2 mag it 
is consistent with being an early B dwarf. 

We conclude from the CMDs and distance estimates. 



that the RSG, the WR star, and the BSGs are all part of 
the same stellar cluster. The average spectrophotometric 
distance of 3.7 ± 1.7 kpc is consistent with the kinematic 
distance 4.7 ± 0.4 kpc within uncertainties. We assume 
the kinematic distance. 

5. THE X-RAY ASSOCIATIONS AND ENERGETICS 

The 2MASS CMD (Fig. ^ shows two distinct popula- 
tions of Chandra sources. One is associated with a blue 
foreground main sequence population, while the other is 
associated with redder and brighter stars that have col- 
ors and magnitudes consistent with a cluster member- 
ship. X-ray sources with a bright 2JVIASS counterpart 
{Ks< 9.6 mag), such as candidate massive members of 
the stellar cluster, are more spatially concentrated; ex- 
cept for one, they are within 3.'5, while the other X-ray 
detections are within 12' from the cluster center. The 
nine candidate massive members with an X-ray associa- 
tion ajeJncluded^nJTable 1. Star #4, the WR star (#8 
in iHadfield et al."2007i ), is associa ted with the Chan- 
dra source #24 ( Hel fand _ et al.ll2007T) . c oinciding with the 
XMJVI source #2 of iFunk et al.l (|2007D . 

IMassive stars can emit X-rays. Single OB stars with 
shocked stellar winds can emit with a typical X-ray lu- 
minosity Lx of 10^1 "33 erg s'^ (|Pollocklll987l ). Shocks 
between the colhding winds of OB-j-OB or OB-I-WR bi- 
narie s can generate a Lx of iQ^^-ss erg s~^ (jClark et al.l 
[2008h . 

To estimate the X-ray fluxes of the nine bright in- 
frared stars associated with an X-ray source, we used 
an interstellar extinction of Ay= 9.1 mag, which was 
estimated from the CMDs, as well as the relationship 
between extinction and total hydrogen column density 
N{H)[cm-^] = Av[mag] x 1.8 x 10^^ Count rates 
were estimated from the Chandra/ACI S-I counts (in 
the b and 2-10 kev) that were given by iHelfand et al.l 
(|2007f ). We converted them into X-ray unabsorbed fluxes 
by using the PIMMS v3.9d, an ACTS count estimator, 
with a power law ni odel and a photon index of 1.5 
(jTownslev et al.l 120061 ). Fluxes were converted into X- 
ray luminosities Lx by assuming a distance of 4.7 kpc. 
Six of the sources have Lx= 2 — 10 x 10^^ erg s~^. The 
three remaining sources (#4, #5, and #10) are an order 
of magnitude brighter (2 — 8 x 10"^^ erg s~^); the bright- 
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est, #4, is the WR star (jHadfield et al.ll2007f) . The ra- 
tio between the X-ray and bolometric luminosities was 
found to be 3.8 x 10~^ for the WR star, but about 10 
times fainter (0.4 x 10^^) for #2 and #3, as expected 
for stars with spectral type later than Bl (jCohenl [19961 
IWaldron fc Cassinellil 120071 ) . We conclude from the X- 
ray energetics that these bright infrared stars are likely 
to be massiv e stars. The WR is most lik ely a colliding 
wind binary (jPortegies Zwart et al.ll200^ . 

6. AGE AND MASS 

The large variety of evolved objects - 1 WR, 1 RSG, 
2 BSGs, and several X-ray emitters - allows us to 
constrain the age and mass of the stellar cluster by 
assuming coevality. The non-rotating Geneva mod- 
els with solar abundance predict the onset of RSG 
stars at an approximate age greater than 6 — 7 Myr 
(jLeieune fc Schaerei1l2001h . In a population of 6 Myr the 
most massive st ars have initial masses Mjn itiai of approx- 
imately 29 Mq (jLejeune fc Schaereill200ll ). In contrast, 
WR stars are present in populations younger than 7.9 
Myr (Minitiai> 20 M©). In particular, single WR stars 
have initial masses greater than 26 — 30 Mq , whereas bi- 
nary WR stars have initia l masses greater than 20 — 25 
Mq (jEldridge et al.ll2008[ ). The X-ray emission associ- 
ated with our WR star suggests a binary system. We 
conclude that the cluster is 6 — 8 Myr old since this age 
allows for the coexistence of both WR and RSG stars. 
Models with non-zero rotation increase both limits by 
-20% (jMevnet fc Maedeil[2005h . 

Assuming that the other eight X-ray emitters asso- 
ciated with the cluster, other than the WR star, are 
BSGs with masses larger than 20 Mq, and by assuming 
a Salpeter IMF down to 1.0 M©, we derive a total initial 
cluster mass of 2000 Mq. If we add as potential BSGs 24 
other stars within 3f5 of the center and with colors and 
magnitudes similar to the X-ray sources {Ks< 9.6 mag, 
J—Ks< 3 mag), the number of BSGs increases to 33, and 
the estimated initial mass of the cluster to 6500 Mq . 

7. SUPERNOVAE REMNANTS 

The cluster core is about 4f5 away from HESS 
J1813— 178. Several studies have recently been carried 
out in order to unveil the nature of the HESS 7-ray 
source. As a result, numerous radio and high-energy 
sources have been detected (see Fig. [T]). Two non- 
thermal radio shells, G12.82-0.02 and G12.72-0.00 
were identified on the 90 cm VLA survey (jBrogan et al.l 
l2006f) . a few arcminutes one from another (see Fig. [T]). 
Very little is known about the SNR G12. 72-0.00. The 
other SN R was associated wit h the HESS J1813-178 
source by iHelfand et al.l (|2005D . A pulsar wind nebula 
(PWN) was detected by Chandra within this SN radio 



shell less than 1' from the maximum probability cen- 
troid of the HESS source (jHelfand et al.ll2007l : lFunk et all 
|2007() . By analyzing the X-ray flux, the authors con- 
cluded that the SN G12.82-0.02, the HESS J1813-178 
source, and the PWN lie at or just beyond 4 kpc, and 
might be associated with the star forming region W33. 
Similar conclusions were also obtained with the XMM 
data bv lFunk et al.l (|2007D . No stellar counterpart to the 
PWN was found in the 2MASS or Spitzer/GLIM PSE 
images bv' Helfand et all (|2007t ). iDean fc Hilll (|2008D as- 
sumed a distance of 4.5 kpc and an age of 300 yr, and 
found the putative pulsar to have an extremely high mag- 
netic field (B = 1.28 x lO^"* G_). iFunk et all (l2007l ) de- 
tected six other X MM sources in additio n to the PWN, 
AXJ1813-178, and IHelfand et al.l ()2007f ) detected a to- 
tal of 75 Chandra sources in the region surrounding the 
twoSNRs (Fig.[T]). 

The new stellar cluster is coincident (to within 1.'6) 
with the radio shell of SN G12.72-0.00, which suggests 
its association with the supernovae progenitor. The kine- 
matic cluster distance of 4.7 kpc is consistent with the 
distance to the high-energy source HESS J1 813— 178, as 
inferre d fro m the hydrogen column density bv lFunk et al.l 
(|2007f) and IHelfand et all (|2007t ). Massive members of 
this cluster were most likely the progenitors of the two 
supernovae and of the pulsar associated with HESS 
J1813— 178. The progenitors of these objects had likely 
an initial mass similar to that of the RSG and WR stars 
(20-30 Mq). 
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